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ABSTRACT 

Using deep high-resolution multi-band images taken with the Very Large Tele- 
scope and the Hubble Space Telescope, we discovered a new anomalous sequence 
in the Color Magnitude Diagram of uj Cen. This feature appears as a narrow, 
well-defined Sub Giant Branch (SGB-a), which merges into the Main Sequence 
of the dominant cluster population at a magnitude significantly fainter than the 
cluster Turn-Off (TO). 

The simplest hypothesis assumes that the new feature is the extension of the 
anomalous Red Giant Branch (RGB-a) metal-rich population discovered by Lee 
et al. (1999) and Pancino et al. (2000). However, under this assumption the 
interpretation of the SGB-a does not easily fit into the context of a self-enrichment 
scenario within uj Cen. In fact, its TO magnitude, shape and extension are 
not compatible with a young, metal-rich population, as required by the self- 
enrichment process. The TO level of the SGB-a suggests indeed an age as old 
as the main cluster population, further supporting the extra-cluster origin of 
the most metal rich stars, as suggested by Ferraro, Bellazzini & Pancino (2002). 
Only accurate measurements of radial velocities and metal abundances for a 
representative sample of stars will firmly establish whether or not the SGB-a is 
actually related to the RGB-a and will finally shed light on the origin of the 
metal rich population of uj Cen. 

Subject headings: Globular clusters: individual (cj Cen); stars: evolution - stars: 
Population II 
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1. Introduction 

The observational facts collected over the last 40 years have demonstrated that the most 
massive (M ~ 5 x 10^ ^o) star cluster in the Halo of the Galaxy, u Centauri, also shows 
the most surprising properties ever observed in a globular cluster. It is the only known 
Galactic globular with a clear metaUicity spread (Norris & Da Gosta, 1995; Smith et al., 
2000). Even more interesting, the abundance variations appear to be correlated with some 
pecuhar structural and kinematical properties (Norris et al. 1997; Pancino et al. 2003). 

Recent wide field photometric studies (Lee et al. 1999; Pancino et al. 2000) have 
identified a previously unknown anomalous RGB sequence (RGB-a), clearly separated from 
the bulk of the RGB population. The RGB-a photometric properties suggest a very high 
metal content, confirmed by spectroscopic follow-up ([Fe/H]~ —0.6; Pancino et al. 2002; 
Orighaet al. 2003). 

In a self-enrichment scenario, the RGB-a would belong to the last burst of star formation 
and should thus be younger than other populations in uj Gen. However, using the photometry 
by Pancino et al. (2000) and the proper motions of van Leeuwen et al (2000), Ferraro, 
Bellazzini & Pancino (2002) showed that the RGB-a stars have a coherent bulk motion with 

respect to the other cluster stars, suggesting that the RGB-a could be a captured component, 
now trapped in the potential well of the cluster. This evidence could eventually fit into a 
different scenario where the cluster is the remnant (i.e., the nucleus) of a disrupting dwarf 
galaxy (see Bekki & Freeman, 2003), not dissimilar from the Sagittarius case (Majewski et 
al. 2000). 

In this letter, we present the discovery of an additional anomalous sequence in the Color 
Magnitude Diagram (CMD) of uj Gen: a narrow Sub Giant Branch (SGB-a) which merges 



^ Based on observations collected at the Very Large Telescope of the European Southern Observatory, 
Cerro Paranal, Chile, within the observing programme 68.D-0332. Also based on observations with the 
NASA/ESA HST, obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., 
under NASA contract NAS5-26555. 

^Dipartimento di Astronomia Universita di Bologna, via Ranzani 1, 1-40127 Bologna, Italy, fer- 
raro@bo.astro.it 

•^INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 1-40127 Bologna, Italy, anto- 
nio@omega.bo.astro.it, pancino@bo.astro.it, bellazzini@bo.astro.it, origlia@bo.astro.it 

*INAF-Osservatorio Astronomico di Teramo, 1-64100 Teramo, Italy, straniero@te.astro.it 

^Department of Physics and Astronomy, San Francisco State University, 1600 HoUoway Avenue, San 
Francisco, CA 94132, USA, cool@sfsu.edu 



-3- 



into the Main Sequence (MS) of the cluster with a Turn-Off (TO) significantly fainter than 
the dominant cluster population. This new feature adds another piece to the complex cu Cen 
puzzle. 

2. Observations and data analysis 

The photometric data consist of a set of high resolution images obtained with FORSl 
at the VLT. All images were obtained (in service mode) between April and October 2002, 
under very good seeing conditions (FWHM= 0'.'4 — 0'.'9), using the high-resolution FORSl 
imaging mode (O.l'.'/pix) that gives a field of view of 3^4 x 3'4. We covered a global ~ 9f x 9f 
area with a mosaic of 3 x 3 pointings around the cluster center. Each pointing consists of 
four deep B, V and 1 exposures, typically of 300, 45 and 25 s, complemented by a number 
of short exposures (typically a few second long) in order to properly measure the brightest 
stars. 

Similar observations (i.e similar strategy and spatial coverage) were obtained on June 
2002 with the Advanced Camera for Survey (ACS) on board HST through three filters 
[namely F625W (i?625), F435W (-B435), F658N (Ha)] within a programme devoted to the 
search for Cataclysmic Variables in u; Cen. We retrieved three 340 s exposures in the S435 
and Rq25 filters for each of the nine pointings from the ESO/ST-ECF Science Archive. Two 
additional short exposures of 8 and 12 s, in the i2625 and ^435 bands, respectively, were also 
retrieved in order to properly measure the brightest stars. 

Both the FORSl and ACS datasets have been reduced with DAOPHOT II (Stetson 
1987, 1994) using the same procedure. To optimize the search for faint objects, the source 
detection has been performed on the combined median image, obtained from the deepest 
exposures of each pointing. The list of detected objects in the combined image has then 
been fed to the PSF-fitting routines of DAOPHOT and ALLSTAR that were run on each 
exposure, separately. Each ACS pointing was corrected for geometric distortions using the 
prescriptions by Hack & Cox (2001). The resulting instrumental magnitudes have been 
converted into a common (instrumental) photometric system and averaged together. The 
final FORSl and ACS instrumental catalogs have then been photometrically calibrated: the 
former, by using a large sample of stars in common with the WFI catalog by Pancino et al. 
(2000) (see also SoUima et al. 2004 for more details), while the latter has been transformed 
into the VEGA-MAG system by using preliminary zero-points, kindly provided by G. De 
Marchi (private communication). 
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3. Results 

Figure 1 shows the {V,B — V) CMD obtained from the ground-based FORSl VLT 
observations. Photometric errors have been calculated for each star from the r.m.s. scatter 
of repeated exposures. In order to minimize the probability of photometric contamination 
by nearby stars, only the 14,000 stars with the lowest photometric errors (cry ~ as <0.01 
mag, see SoUima et al. 2004, in preparation) have been plotted. The complex structure of 
the RGB shows beautifully, but the most striking feature in this CMD is the additional, 
narrow SGB sequence, which sticks out of the main SGB-TO region of the cluster. 

Fig. 2 displays the (-^625, -B435 — -^625) CMD obtained from the HST/ACS sample: more 
than 400000 stars with the lowest photometric errors are plotted. All the evolutionary 
sequences are extremely well defined, apart from the brightest giants that are saturated even 
in the shortest exposures (not plotted in Fig. 2). The SGB-a is better defined than in the 
FORSl CMD, and the complex SGB/TO structures are more clearly distinguishable. In this 
CMD, the SGB-a appears as the direct continuation of the RGB-a, although the presence 
of a number of stars located in between the main SGB and the SGB-a (at Rq25 = 17 — 18) 
offers room to other possible interpretations. The SGB-a sequence merges into the bulk of 
the cluster MS at a significantly fainter magnitude {AV ~ Ai?625 ~ 0.6 mag) with respect 
to the dominant population. SGB-a stars in common to both the FORSl and ACS catalogs 
are marked in Figure 1. 

Both the FORSl and ACS CMDs show that the MS-TO of the SGB-a population 
{Rjg{SGB - a) = 18.40 ± 0.15) is significantly fainter than the MS-TO of the dominant 
cluster population {Rq^{MP) — 17.80 ± 0.15). A hint of the presence of such a feature was 
already noted by Anderson (2002, see his Fig. 1), in the core of the cluster on the basis of an 
instrumental (C/336, f/336 — -^675) CMD obtained from WFPC2/HST observations. However, 
the high quality CMDs presented here provide the first clear-cut definition of this sequence 
in the CMD of uj Gen. 

The TO and SGB regions of u; Gen show the same level of complexity already observed 
at the RGB level: (i) a well defined and populated main SGB sequence, highly homogeneous 
both in terms of age and metallicity. (ii)) a more continuous distribution of stars toward 
the red, tracing a less homogeneous stellar component, with different metallicities and/or 
ages, (in) a well defined narrow additional SGB, redder and fainter than the main SGB. 
The sharpness of this sequence strongly suggests a high degree of homogeneity in the age 
and metallicity of its stars. 
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4. Discussion 

In the simplest scenario, the SGB-a is the extension of the RGB-a. Following this 
hypothesis, we assume that the anomalous branches are constituted by the same population 
of metal-rich stars. 

We then performed a detailed comparison with theoretical models by using a set of the- 
oretical isochrones calculated adopting the most up-to-date input physics (Straniero, Chieffi 
& Limongi 1998). In particular, the equation of state includes the electrostatic correction 
(see Prada Moroni & Straniero 2002, for a detailed description) and microscopic diffusion 
(gravitational settling and thermal diffusion). The theoretical isochrones have been trans- 
formed into the observational planes by means of the synthesis code described in Origlia & 
Leitherer (2000), using the model atmospheres by Bessell, CasteUi & Plez (1998). For the 
ACS data set the filter responses and camera throughputs, kindly provided by the ACS User 
Support Team, have been used. 

In order to perform a meaningful comparison between theoretical isochrones and ob- 
served sequences, we need to adopt an average metallicity for each population. As discussed 
in previous papers (see Ferraro et al. 1999, 2000 and references therein) an accurate com- 
parison requires the use of the global metallicity, [M/H], obtained by combining the classical 
[Fe/H] abundance with the contribution of the a— elements. Detailed abundances for the 
different stellar populations in a; Cen have been published by several authors (see, e.g., Nor- 
ris & Da Costa 1995; Smith et al. 2000; Vanture et al. 2002), including the recent studies 
on the RGB-a stars published by our group (Pancino et al. 2002, 2003; Origlia et al. 2003). 
These papers show that, while the metal poor populations at [Fe/H]< —1.0 turn out to be 
a-enhanced ([a/Fe]^ 0.3 — 0.4 dex) compared to the solar values, the most metal rich one 
is only marginally a-enhanced [Q;/Fe]< 0.1 — 0.2 dex. 

Using the above values and the prescription by Salaris, Chieffi & Straniero (1993, see 
their eq. 3), we adopted a global metaUicity [M/H]~ —1.5 {ox Z — 0.0006) for the metal- 
poor, dominant population and [M/H]~ —0.6 (or Z — 0.005) for the metal- rich population. 
The mean interstellar extinction coefficients listed in Table 2 by Savage & Mathis (1979) 
have been adopted. 

By best-fitting the Z — 0.0006 isochrone to the metal-poor component, we obtained 
(m — M)o = 13.70 and E{B — V^ — 0.11 and an age of 15Gyr. These values are in good 
agreement with the most recent determinations by Thompson et al. (2001) and Lub (2002) 
((m — M)o = 13.65 ±0.11 and E(B — V^ = 0.11 ±0.02) and the most recent age determination 
by Gratton at al. (2003). By using these values, the position of the anomalous SGB has 
been compared with a set oi Z = 0.005 isochrones in the 12-18 Gyr age range. It is worth of 
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mentioning that, since we are performing a strictly differential analysis between the SGB- 
a and the dominant cluster population, small uncertainties in the reddening and distance 
modulus would not affect the overall results. 

As can be appreciated in Fig. 3, while the metal-poor population is reasonably well re- 
produced, the TO level of the anomalous population is definitely fainter than what predicted 
by a metal-rich isochrone significantly younger {t < lAGyr) than the dominant population. 

A detailed analysis of Fig. 3 shows that indeed the observed morphology of the anoma- 
lous SGB is not correctly reproduced by any Z — 0.005 isochrone, regardless of age. In fact, 
though the MS-TO level of the anomalous population appears consistent with an age of 17 
Gyr, (i.e. 6t ^ 2 Gyr older than the metal-poor component), the TO color is significantly 
bluer than the corresponding isochrone. Moreover, the morphology, extension and position 
of the lower RGB-a and SGB-a are not correctly reproduced, since 1) the SGB-a appears 
significantly less extended in color and much steeper than what predicted by a metal rich 
isochrone (somewhat suggestive of a lower metal content); 2) the position of the base of 
the RGB-a appears significantly bluer than what predicted by the isochrone.^ The same 
problems have been noted by fitting the SGB-a in the (V,B-V) plane. Note that the same 
isochrone set nicely fits the SGB/TO region of well-known clusters of similar metaUicity (as 
47 Tuc). 

Other two parameters, namely the He and CNO abundances can affect the location and 
shape of the SGB-TO region in the CMD. In order to evaluate the effect of an enhanced 
He abundance over the shape of the SGB, we computed a set of suitable models at Y=0.28 
{SY — 0.05 with respect to the standard value adopted above). From this models we noted 
that an increase of the He abundance does not significantly affect the TO level while it 
decreases the SGB extension by moving blueward the location of the RGB by few hundredths 
of magnitude in the B — R color. Hence an increase of the He abundance cannot be invoked 
to explain the faintness of the TO of the anomalous population. 

Renzini (1977, see also the discussion in Salaris, Chieffi & Straniero 1993) showed that 
the TO luminosity and temperature primarily depend on the CNO and Ne abundances, 
while the RGB colors mainly depend on the Fe, Si and Mg ones. Thus an ad hoc increase 
of Zc,N,o,Ne would (i) movc the TO color towards the red while leaving the RGB location 
unchanged, thus reducing the SGB extension and (ii) move the TO towards fainter magni- 
tudes (see Fig. 22 of Salaris, Chieffi & Straniero 1993). This goes exactly in the direction 
of moving the young {^12 Gyr) isochrone towards the observed sequence. However, Salaris, 



^However, when cvahiating the importance of this discrepancy one needs to keep in mind that the position 
of the RGB in the theoretical models is highly uncertain due to the mixing-length convection approximation. 



- 7- 



Chieffi & Straniero (1993) quantified the effect of O (tlie most abundant element among 
tfie CNO group and Ne) overabundance over the TO luminosity: an O overabundance by a 
factor of 4 would increase the TO level by SV ~ 0.07 mag. An unlikely O overabundance 
of ~ 20 would be necessary to move the 12 Gyr isochrone to the level of the observed TO. 
Note that spectroscopic [O /Fe] measurements of a few RGB-a stars (see Fig. 5 in Origlia et 
al. 2003) indicate only marginal (if any) O overabundance. 

In summary, it seems impossible to reconcile the TO level of a young (^ 12 Gyr) 
metal rich isochrone with the observed feature. Accurate measurements of C,N,Ne and O 
abundances for a significant number of SGB-a stars are urgently needed to check for possible 
anomalies, which could play some role in characterizing the morphology and the position of 
the SGB-a. Moreover, additional modeling with a detailed fine-tuning among [M/H], CNO 
and Ne abundances could indeed improve the agreement between old (14-16 Gyr) isochrones 
and the observed SGB-a morphology. 

A number of papers (Lee et al. 1999, Hughes & Wallerstein 2000, Hilker & Richtler 
2000, Rey et al. 2003) argue, from indirect evidences, that the age of the most metal rich 
stellar population in a;Gen is significantly {At ~ 4 Gyr) younger that the most metal poor 
one. Under the hypothesis that the SGB-a is related to the RGB-a, we demonstrated from 
the direct detection of the TO, that this anomalous, metal rich population cannot be younger 
than the most metal poor one. This evidence excludes the possibility that such a metal rich 
population be the latest product of the self-enrichment process in u Ccn (as suggested by 
Rey et al. 2003 and references therein), while it would support its extra-cluster origin (see 
Ferraro, Bellazzini, Pancino 2002). 

In the scenario proposed by Bekki & Freeman (2003), uj Cen is the complex stellar relic 
of a nucleated dwarf galaxy that merged in a remote epoch with the Galaxy (see also Dinescu 
et al. 1999; Majewski et al. 2000). Populations with different metallicities would be the 
result of subsequent radial gas inflows into the main body of the parent dwarf galaxy. This 
scenario can reasonably account for the genuine cu Cen sub-populations at [Fe/H]< —1.0. 
Conversely, the anomalous, metal rich population could be a small stellar system (a globular 
cluster?) accreted by the main body of cu Cen during the disrupting interaction with the 
Galaxy. The assumption that oo Cen could host a satellite globular cluster is not fully ad 
hoc, since both the other two known disrupting systems in the Galaxy (the Sagittarius dwarf 
spheroidal - Ibata et al. 1994 and the Canis Major dwarf galaxy - Ibata et al. 2003) have 
their own globular cluster system. 

A number of alternative hypothesis about the nature of the SGB-a become possible 
if this population is not related to the RGB-a. The SGB-a observed morphology could 
be indeed reproduced by old, low-metallicity isochrones, but this would require the use of a 
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different distance and reddening with respect to main cluster population. Such an hypothesis 
would imply that the anomalous SGB trace a metal-poor stellar system, in the background 
of the main body of uj Cen. Presently we have no possibility to discriminate between the 
various scenarios. Only direct measurements of radial velocities and metal abundances of a 
significant sample of SGB-a stars will allow us to answer the enigma posed by the discovery 
of this new feature. 
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also thank the Referee, Russell Cannon, for his comments and suggestions. This research 
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Fig. 1. — (V, B — V) CMD for stars identified in tlie nine fields observed with FORSl at the 
ESO/VLT. About 14,000 stars with the lowest photometric error have been plotted. The 
SGB-a population, merging into the cluster MS, is clearly visible. Stars marked with large 
open circles are identified as SGB-a stars in both the FORSl and ACS catalogs. 
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Fig. 2 — (-^625,-8435 - -R625) CMD for more than 400,000 stars identified in the nine ACS 
fields. The SGB-a is clearly visible, along with the complex sub-structures of the TO-SGB 
region of cu Cen. 
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Fig. 3. — Zoomed (i?625,-B435 — -^625) CMD of the TO region. Two isochrone sets are 
overplotted. The metal-poor dominant population is well reproduced hy a, Z — 0.0006 
isochrone of 15 Gyr (heavy sohd line). A set oi Z — 0.005 isochrones in the 12-18 Gyr age 
range is also shown. 



